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Abstract 

The method of spectral decimation is applied to an infinite collection 
of self-similar fractals. The sets considered belong to the class of nested 
fractals, and are thus very symmetric. An explicit construction is given 
to obtain formulas for the eigenvalues of the Laplace operator acting on 
these fractals. 

In 1989, J. Kigami [5] gave an analytic definition of a Laplace operator act- 
ing on the Sierpinski Gasket; a few years later, this definition was extended to 
include Laplacians on a large class of self-similar fractal sets [6], known as post 
critically finite sets (p.c.f. sets). The method of spectral decimation introduced 
by Fukushima and Shima in the 1990's provides a way to evaluate the eigen- 
values of Kigami's Laplacian. In general terms, this method consists in finding 
the eigenvalues of the self-similar fractal set by taking limits of eigenvalues of 
discrete Laplacians that act on some graphs that approximate the fractal. The 
spectral decimation method was applied in [1] to the Sierpinski Gasket, in order 
to give an explicit construction that allowed to obtain the set of eigenvalues. In 
[TT] it is shown that it is possible to apply the method to a large collection of 
p.c.f. sets, which includes the nested fractals defined in 0. In addition to the 
Sierpinski Gasket, several specific cases have been treated in the literature (e.g. 

EH HI El El). 

In the present work, we develope in an explicit way the spectral decimation 
method for an infinite collection of self-similar sets, parametrized by n > 2 (a 
positive integer) . The definition of these sets is given in Definition [l] The cases 
n = 2, 3, corresponding respectively to the unit interval and to the Sierpinski 
Gasket, have been presented with thorough detail in [TU]. Our presentation 
follows this reference to some extent. However, some technical difficulties arise 
for n > 4. This is mainly due to the fact that -even though the fractals con- 
sidered here are still very symmetric - the graphs approximating the fractal are 
not as homogeneous as the ones for the Sierpinski Gasket. For instance, if we 



Figure 1: The approximating graph Tx to the left (for P 4 ), and 1^ without the 
boundary to the right. 

consider the graph obtained by by taking away the boundary points from Fi 
(see Definition [2] and Figure [l}, then it will be a complete graph only for n < 3; 
a consequence of this, is the appereance of sets of vertices that have to be dealt 
with (and which we denote by G r , s ) that are not present for n < 3. 

In Section [l] we present general facts about self-similar sets, for the sake of 
completeness and in order to establish notation; at the end of the section we 
introduce the sets P n of particular interest in this work. In Section [5] we define 
the graphs that approximate the self-similar sets P„ and fix more notation. 
Our main result is developed in Section [3] (Theorem [T]); it is shown that the 
eigenvalues and eigenfunctions of the discrete Laplacians of the approximating 
graphs can be obtained recursively. Finally, in Section |4j it is shown that the 
eigenvalues of the Laplace operator in P„ can be recovered by taking limits of 
the discrete Laplacians; in order to do this, we solve the renormalization problem 
for this case (see Theorem [2]). 

1 Notation and Preliminaries 

We denote by S„ the shift-space with n symbols. In this work we will always 
consider these n symbols to be the numbers 0, 1, . . . , n — 1. S n is a compact 
space (see e.g. [7]) when equipped with the metric 

<5(a aia 2 a 3 . . . ; b a b 1 b 2 b 3 ...) = — 
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where 

k = min{j > | a,j ^ bj}. 

We will use the dot notation a, meaning that the symbol a repeats to infinity. 

Let x = X0X1X2 ... an element of S„, and a G {0, . . . ,n — 1}. We denote by 
T a the shift-operator given by the contraction 

T a (x) = ax Q Xix 2 ■ . . 

The space S n is a self-similar set, equal to n smaller copies of itself, with 
{T , . . . , T n _i} the corresponding contractions. Even more so, it can be proved 
(Theorem 1.2.3. in 7 ) that every self-similar set equals a quotient space S„/ ~ 
for an equivalence relation ~. 

For a = (aoai . . . a m _i) a word (in the alphabet {0, . . . , n — 1}) of length m, 
we denote by T a the shift-operator given by 

T a (x) = a ■ ■ ■ a m -iXoX-L%2 ■ ■ ■ 

and call it an m- contraction; the sets of the form T a (S„) are known as cells of 
level m. We note that, for each choice of m, S n is the union of the n m cells of 
level m. 

The self-similar fractals that are considered in this work arc defined next. 

Definition 1. For n G N define P„ as the quotient space S„/ ~, with the 
equivalence relation given by 

aoaia2 . . • Ofc&c ~ aoaia2 . . . a^cb 

for any choice of symbols a,j, b and c. 

Pi is a trivial space with only one element, Pi is homeomorphic to a compact 
interval in K, and P3 is the Sierpinski Gasket. 

2 Graph approximations of Self-Similar Sets 

In this and the next sections, we consider one of the self-similar sets P n defined 
above, for an arbitrary but fixed value of n > 2. 

Let Vq be the set of points in P„ that have the form k with fc = l,...,n — 1. 
We call Vo the boundary of P„. Likewise, for m € N let V m be the subset of 
P„ of points of the form oq . . . a m -xk. In other words, x G V m if and only if it 
belongs to the image of Vo under some m-contraction. 

Next, we define the graphs that will approximate P„. 

Definition 2. Denote by Tq the complete graph of n vertices, with Vq its set 
of vertices. For m £ N, let F TO be the graph with set of vertices V m and edge 
relation established by requiring x to be connected with y if and only if there 
exists an m-contraction T a such that both points x and y are in T a (Vb). 
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We can see that an equivalent formulation is that two vertices x and y share 
an edge in T m only when their first m symbols coincide. It is worth noting that 
even though Vq C V\ C Va • • • , the edge relation is never preserved; this follows 
from the fact that if x 7^ y are connected in T m , then their (m + l)-th symbols 
cannot be equal, so that they will not be connected in r m+1 . 

For each m € N let A m be the graph-Laplacian on T m . We consider the 
Laplacian as acting on a space with boundary. More precisely, for a real-valued 
function u defined on V m and x in V m \ Vq: 

A m u(x) = ^2(u(x) - u(y)), 

with the sum over all vertices y that share an edge with x\ the boundary values 
remain unchanged. Also u is an eigenfunction of Ao with eigenvalue A, if 

A m u(x) = Xu(x), Vie € V m \ Vq. 

We denote by E m (-, •) the associated quadratic form (known as the energy prod- 
uct of the graph): 

E m (u,v) = (A m u,v) = y^(u(x) - u(y))(v(x) - v(y)) 

for u and v real-valued functions defined on V m , and the sum being taken over 
the pairs of vertices (x,y) that are connected to each other. Also, we use the 
abbreviation E(u) = E(u,u). 

3 Spectral Decimation 

Let m > 1, and suppose u is an eigenfunction of A m _i, with eigenvalue A m _i. 
We will show that it is always possible to extend this function to the domain 
V m so that it will be an eigenfunction of A m (with not the same eigenvalue). 
In order to do this, we will derive necessary conditions for the extension to be 
an eigenfunction; in the process, it will become clear that those conditions are 
also sufficient. 

Suppose that u is an eigenfunction of A m with eigenvalue A m ; we aim to 
write the values of u m in V m \ V m -\ in terms of its values in V m -\. Without 
loss of generality, we can restrict ourselves to the set V m n T a (P n ) for a fixed 
(m — l)-contraction T a ; this is because the vertices of T m that belong to the set 
(V m \ V m -i) n T a (P„) are not connected to any vertices outside the cell T a (P„). 
Denote the elements of this set by 

Xf, jC = abc. b, c = 0, . . . n — 1. (1) 

It is clear that Xb.c = %c.b, and also that Xb, c € V m -i if and only if b = c. 
This is shown in Figure [2j 
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Figure 2: A cell of level m — 1 of a graph T m , approximating P4. 

For each point x r , s £ V m n T a (P„) define the sets of vertices 

^r, s = {Xi,j {r, s} n {», j} 7^ 0} 

G r , s = {a^ I i ^ j, {r, s} n {t, j} = 0} 

This is, F r!S is the set of vertices (not in r m _i) that are connected to the vertex 
x r>s in T m , and G r ^ s is the set of vertices (not in r m _i, either) that are not 
connected to it. 

For every r/swe have 

(2(n - 1) - \ m )u{x r ^ s ) = u(x rtr ) + u(x s , s ) + ^ u(xij). (2) 

F r ,s 

Adding this up over all the possible values of r and s, and rearranging terms 
yields 

n-1 

(2 - A m ) ^ u ( x r,s) = (n - 1) ^ u(xjj), 
which for any fixed a ^ b can also be written in the form 




This, together with ([2]) allows us to express the sum of the values in G a ,b in 
terms of u(x a ^) and the values at points in V m -\\ namely, provided A m 7^ 2, we 
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have that 



(n — 1) n— 1 

Z — Am 

G a , b m j=0 

Next, we will take the sum of the same terms, but only over the Xi.j inside 
the set F a ^ b for fixed values a =^ b. In order to do this, we note that F ab contains 
two complete graphs with n — 1 vertices (one with the points Xi^ a and one with 
the points cc^b), and these complete graphs are pairwise connected to each other 
(xi.b ~ Xi.a). From this, it is clear that each Xij € F ab is connected to other 
n— 1 vertices in F a<b . Also, each Xij € G a ,b is connected to exactly four vertices 
in F a: b- For the vertices in V m -\ we note that Xjj is connected to n — 2 vertices 
in Xij € F a b if j = a,b and to only two vertices otherwise. 

From the preceeding discussion it follows the equality 



(n - A m ) ^ u ( x i,o) = 4 X! u ( Xi d) + ( n - 2)(u(x a , a ) + u(xb,b)) 

Fa,b G a b 

+ 2(n - 2)u{x afi ) + 2 u(x jd ) (4) 

Consider the expresion given by ^ for {a, b} = {r, s}, multiply it by n — A m , 
and substitute equality Q into it; this gives after arranging terms 



(X 2 m - (3n - 2)A m + 2(n 2 - 2n + 2))u(x a ,b) = 4 ^ u ( x i,i) 

G a , b 

+2 2J u(xj,j) + (2(n - l)-\ m )(u(x ata ) + u(x b ,b))- 

j¥=a,b 

We want to get rid of the terms corresponding to G a>b , so we replace it by pH). 
After straightforward computations, we can see that for A m 7^ 2 

(A^ - (3n + 2)A m + 2n{n + 2))u(x«, 6 ) = 

(A^ - 2{n + 2) + 8n) (u(a; a , a ) + u(x bib )) + 2{2n - A m ) "(^jj) 

2 — A m 

The quadratic equation for A m in the left hand side has roots n + 2 and 2n. 
The one in the right hand side has roots 4 and 2n. This gives us the following 
expression for u(x a ,b) m terms of the values of u in V m —\: 

_ (4-A m )(u(y) + u(i„))+2[; ]/rs B(i ]J ) 

(2-A m )((n + 2)-A m ) (5j 

valid for any eigenvalue A m 7^ 2, n + 2, 2n. 
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For A m = this reduces to 

2 ^ , 

u(x r ,s) = rn( u ( x r,r) + u(x s . s )) H — - } u(x jt j). (6) 

n + 2 n + 2 ^— ' 

It is clear from the construction that if u{x a ^) is defined by ^ and A m is 
given by (10 1, then we have that 

A m it(a; Qj 6) = X m u(x a>b ). (a^b). 

It remains to verify that this is valid as well in V m -i- Of course, this cannot 
be true for arbitrary values of A m , but only at most for specific values depending 
on A m _i; we will find those values in what follows. 

Take a point in V m -\, say 

x p,p — a P: a = Oq . . . a m -i- 

Suppose that = q is the last symbol in a that is different from p; we can 
assume that such symbol exists, since otherwise x p p would be in the boundary 
Vq. With this, the point x PtP can also be written in the form 



p,p 



aq, a = Oq . . . a k -ipq . 



with the necessary number of g's to make a' a word of length m — 1. Hence, x PiP 
is in exactly two different (m — l)-cells: T a (P„) and T a /(P n ), corresponding to 
each one of its two representations. 

Denote by x' r s the points in T a /(P„) n V m , defined as in ([lj for the points 
in T a (P„) n V m ] in particular x PtP = x' q q (see Figure [3]). The value of u in the 
points x' ra is given by the analogue of equation The vertex x PiP is connected 
in T m to the 2(n — 1) points of the form x p j and x' q j, from which it follows that 
u is an eigenfunction of A m with eigenvalue A m , if and only if ^ holds for all 
x r ,s G V m \ V m -i and the following equality holds for all x p , p € V m -\\ 

(2(n - 1) - A m )tt(afp, p ) = ( u (xp,j) + u { x q,j)) ■ ( 7 ) 

On the other hand, since we know that u is an eigenfunction of A. m _i with 
eigenvalue X m -i, we also have that: 

(2(n - 1) - A m _i)u(a;p i p) = + u ( a; i,i) • ( 8 ) 

Replacing each term in the right hand side of ([7| by its expression given by 
([5| we can see that 

11 1) XmMx p>p) - ( 2 -A m )((n + 2)-A m ) 



7 




Figure 3: Two cells of level m — 1 intersecting in a vertex of T ri 



and using Q this gives 

[2(n - 1)(4 - A m ) + (2n - A m )(2(n - 1) - A m )] 



(2(n - 1) - A rn )(2 - A ro )((n + 2) - A m ) 



Taking m(x P;P ) 7^ 0, and cancelling out, after computations the above equal- 
ity reduces to the quadratic 

A^-(n + 2)A m + A TO _i = 0, (9) 

which in turn gives the following recursive characterization of the eigenvalues: 



Ar 



(n + 2) ± N /(n + 2) 2 -4A m _i 



(10) 



Since this procedure can be reversed, we have proved the following result. 

Theorem 1. Let A m 7^ 2,n + 2,2n, and let A m _i &e given by Suppose 
u is an eigenfunction of A m _i wi£/i eigenvalue \ m ~i- Extend u to V m by 
Then u is an eigenfunction of A m luii/i eigenvalue X m . Conversely, if u is an 
eigenfunction of A m with eigenvalue A TO 7^ 2,rt + 2,2n, then its restriction to 
V m -i must be an eigenfunction of~V m —\ with eigenvalue X m -i for either choice 
in ([TO}. 
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4 The Laplacian on the Self— Similar Fractals 



In order to define the Laplace operator of a p.c.f. fractal by means of graph 
approximations, it is required to solve the so called renormalization problem 
for the fractal (e.g. [TO], Chapter 4); roughly, this consists in normalizing the 
graph energies in T m in order to obtain a self-similar energy in the fractal by 
taking the limit. This can be achieved if the energies are such that they remain 
constant for each harmonic extension from T m to r m+ i. Below, we do this for 
the P„ sets. 

Definition 3. For a given function u with domain V m -i, we call the extension 
of u to V m given by ^ its harmonic extension. 

The next result, gives the explicit solution of the renormalization problem 
the for P„. 

Theorem 2. Let u : V m -i — > K arbitrary, and let v! : V m — > K be its harmonic 
extension. Then 

n 

E m {u') = — — -E m -x(u) 
n + 2 

Proof. Note that the energy at level A; of a given function equals the sum 
of the energies at all the /c'-cells for any k! < k, since different cells share 
no edges. This allows to restrict ourselves to one fixed m — 1-cell both while 
considering E m (u') and E rn -i(u). We use the notation of the previous section 
for the vertices of T m in this cell, and write E for the energy restricted to this 
cell. We can readily see that 

E m ^i{u) = 22(u(x iti ) - u(x 3 j)) 2 . 

n-l 

= (n - 1) u 2 (x iti ) - 2^2u(xi i iu(xi,j)). (11) 

i=0 i^j 

In order to evaluate the energy E m we consider first the edges joining vertices 
in V m ~i with vertices in V m \ V m -x- The edge joining the vertex x a ^ a with the 
vertex x a ,k contributes to the energy by 

(u(x a:a ) ~ u(x a ,k)) 2 = , ^ ^ I nu(x a ^ a ) - 2u{x k ,k) ~ ^ u {^j,j) 

When adding up this over all possible pairs a/lt, each x r ^ r will appear n — 1 
times as the x a , a , another n— 1 times as the Xk,k and (n— l)(n— 2) as one of the 
Xj.j's. Each pair double product 2x r-r x StS will appear twice for {r, s} = {a, fc}, 
2(n— 2) times for {r, s} — {a, j} for some j, also 2{n— 2) times for {r, s} — {k,j} 
for some j, and finally (n — 2)(n — 3) times when both r and s are one of the j's. 
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All this implies that the contribution to the energy from these edges is, after 
simplification: 



^(u{ x a,a) - U(x a}k )f 



(n - l)(?i 2 + n + 2) 
(n + 2) 2 



2/ \ 2(n 2 +n + 2) ^ 

<=o ln + ^ j 



(12) 



On the other hand, the contribution from the edge that joins the vertices x a ^ 
and x a , c (in \ V^-i) equals 



(u(x a ,b) - W(^a,c)) 2 



{u(x b . b ) - u(x c . c )Y 



(n + 2) 2 

Taking the sum over all of the vertices in V m \ V m -\ yields 
^ (u(x a>b ) - u(x a , c )) 2 

- ''' 1>(/i ^M 2 (xi,i)- ^^^]u(a:i,i)tt(a:jj). 



(n + 2) 2 



i=0 



From ( 12 ) and ( 13 1 it follows that the total energy of the cell is 

n-l 



E m (u') = 



n(n — 1) 
n + 2 



u 2 (x i4 ) - 2n(n + 2). 



(13) 



(14) 



This, together with (lit gives 

~ Ti ~ 

E m {u') = — —E m -i(u). 
n + 2 

Taking this result for all the m — 1-cells concludes the proof. 



Definition 4. The energy in P n is given by 
E( u )= km f — 

m— >oo \ n 



E m (u). 



a 



The domain of E(-) being the space of functions such that the energy is 

finite. Write for the subspace of of functions that vanish on the 

boundary. The energy product E(u, v) can be recovered by the polarization iden- 
tity. 



Let n be a self-similar measure in P n , the Laplacian is given by: 
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Definition 5. (Kigami's Laplacian) With /i and A p as above, we say that u 
is in the domain of A^ if there exsists a continuous function f such that 



E(u, v) = - J fvdfx Vw e D ( n) 



In such case, we define A„u = f. 



Aside from the above weak representation, a pointwise formula can be ob- 
tained for A m u, proceeding in exactly in the same way as in [10] (Theorem 
2.2.1). In the case where /_t is the standard measure in P n (i.e. the only Borel 
regular measure such that the measure of every m-cell is equal to n~ m ), the 
pointwise formula is 

71 

Aau(x) = ~ lim (n + 2) m A m u(x). 

2 m— >oo 



This leads to the following: If a sequence {A m } is defined recursively by (10 1 
(assuming that A TO is never equal to n, n + 2 or 2n), and u m is given by relation 
((5| then 

— 77 

A = - lim (ti + 2) m A m 

2 m— foo 

is an eigenvalue of A M with eigenfunction u given by the limit u m — > u. The 
limit above exists provided that the sign in relation ( 10 ) is chosen to be "+" for 
at most a finite number of times. 
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